Multiboson spin- wave theory for Ba 2 CoGe 2 07, a spin-3/2 easy-plane Neel 
antiferromagnet with strong single-ion anisotropy 
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We consider the square-lattice antiferromagnetic Heisenberg Hamiltonian extended with a single- 
ion axial anisotropy term as a minimal model for the multiferroic Ba2CoGe2 07. Developing a 
multiboson spin- wave theory, we investigate the dispersion of the spin excitations in this spin-3/2 
system. As a consequence of a strong single-ion anisotropy, a stretching (longitudinal) spin-mode 
appears in the spectrum. The inelastic neutron scattering spectra of Zheludev et al. [Phys. Rev. 
B 68, 024428 (2003)] are successfully reproduced by the low energy modes in the multiboson spin- 
wave theory, and we anticipate the appearance of the spin stretching modes at «4meV that can be 
identified using the calculated dynamical spin structure factors. We expect the appearance of spin 
stretching modes for any S > 1/2 compound where the single-ion anisotropy is significant. 

PACS numbers: 75.30.Ds 75.30.Gw 75.10.Jm 
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I. INTRODUCTION AND THE MODEL 

The spin excitations in isotropic Heisenberg models 
with long-range Neel order are well understood within the 
framework of the spin-wave theory, even in the extreme 
quantum spin- 1/2 case on the square lattice, as exempli- 
fied by La 2 Cu04.— ^ In the spin- wave approach, the spin 
operators are transcribed using a Holstein-Primakoff bo- 
son, and expanding in 1/S around the classical S —> oo 
ground state, one gets a Hamiltonian quadratic in the 
bosonic operators that describes the spin excitations of 
the ordered state. Interestingly, the spin-waves as de- 
scribed above do not take account of all the possibili- 
ties for larger spins. In S > 1/2 quantum spin mod- 
els higher order tensor interactions (like biquadratic ex- 
change) and anisotropy terms (like axial anisotropy) may 
lead to multipolar or nematic ordering^— When studying 
the dynamical properties of such systems, however, the 
conventional spin-wave approach fails, and one needs to 
introduce generalized bosonic operators related not only 
to the spin but also to higher order spin operators^— 
Similarly, additional bosons are needed for spin systems 
with orbital degeneracy— Interacting spin multiplet sys- 
tems have also been studied by an extended Holstein- 
Primakoff theory in which a boson is introduced for each 
energy level of the multiplet. Based on this approach 
the excitation spectra of Cu 2 Fe 2 Ge 4 0i3 and Cu 2 CdB 2 6 
have been theoretically reproduced and for both materi- 
als longitudinal modes have been reported that are re- 
lated to the deformation of the spin wave-function on 
the magnetic ions, also leading to the reduction of the 
magnetic moment of the spinAi A multiboson approach 
with spin-orbital coupled 27 single ion levels has been ap- 
plied to the case of the Lai.sSro.sCoO^ and La 2 CoO^, 
where the Co ions are in octahedral environment. 

Our work is inspired by the strongly anisotropic spin- 
3/2 multiferroic material, Ba 2 CoGc 2 07. This compound 
has tetragonal symmetry and can be characterized by 



layers of square lattices formed by the magnetic Co 2+ 
ionsJ^— As the neighboring cobalts are positioned in 
differently oriented tetrahedral environments, the unit 
cell contains two of these. Below the transition temper- 
ature Tnoei = 6.7 K, the magnetic moments are antifcr- 
romagnctically aligned in the plane of the Co 2+ ionsJ^ 
Spin excitations have been studied by inelastic neutron 
scattering in Ref. ITil . and the observed dispersions were 
fitted using the conventional spin-wave theory based on 
large exchange anisotropy. Additional spin excitations 
at ps 1 THz e nerg ies were observed in light absorption 
spectra in Ref. [l9|. While these additional modes are be- 
yond the conventional spin- wave theory, they were repro- 
duced in an exact diagonalization study of small clusters 
by Miyahara and Furukawaii using a Heisenberg Hamil- 
tonian extended with a strong single-ion anisotropy of 
the (S z ) 2 form. Recently, we applied the aforementioned 
multiboson spin-wave theory to describe the modes ob- 
served in far-infrared absorption spectra in external mag- 
netic fields In the present paper, we aim to study the 
effect of the strong single-ion anisotropy on the zero-field 
spin-wave spectrum in the momentum space using the 
multiboson spin-wave approach, and to give predictions 
for the dispersion of the higher energy modes in the in- 
elastic neutron scattering spectra. 

Following Ref. [l?], we consider the Hamiltonian 



T~L — J ^ ^ (jSf SJ + Sf Sj^J + J z ^ ^ SfSj + 

A£(5) 2 , (1) 



where J and J z are the exchange couplings, A is the 
strength of the single-ion anisotropy, the summation is 
over the (i, j) nearest neighbor sites, and the z-axis is per- 
pendicular to the square lattice plane. In Ref.[l7l a rather 
strong easy-plane anisotropy, A/ J ss 8 has been sug- 
gested in Ba 2 CoGe 2 07. Furthermore, a Dzyaloshinskii- 



Moriya interaction 0.04J has also been considered, 
which is in fact very small compared to the exchange cou- 
pling and anisotropy term, thus we omit it in the present 
study. 

The paper is structured as follows: In Sec.|TT]we shortly 
discuss the variational phase diagram of Hamiltonian ([T]) 
as a function of easy-plane and exchange anisotropies. 
The variational wave function is then used in Sec. IIIII to 
construct a suitable boson basis and perform the multi- 
boson spin-wave approach. The spin-wave Hamiltonian 
is diagonalized numerically and analytically in the mo- 
mentum space, and the behavior of the modes for differ- 
ent exchange and anisotropy parameters is discussed. In 
Sec.|lV]we calculate the dynamical spin structure factor. 
These results are quantitatively compared to the inelastic 
neutron scattering measurements^ of the Ba2CoGe2C>7 
in Sec. |V]and we draw conclusions in Sec. I VI I Finally, in 
Appendix the A — > and A — > oo cases are discussed in 
detail. 



II. VARIATIONAL GROUND STATE 

Let us start with the determination of the ground state 
phase diagram variationally, assuming a site factorized 
trial wave function |\fr) =Y\ i \^i), where the index i runs 
over the lattice sites. The is a wave function in the 
four dimensional local Hilbert space of the 5 = 3/2 spin 
on site The variational phase diagram of the Hamil- 
tonian (TTJ has been discussed previously in Ref.[2l|: in ac- 
cordance with experimental findings^ a two-sublatticc 
antiferromagnetic order is realized for the relevant pa- 
rameters, with j^j) = \^a) if site i is on A sublatticc 
and \^i) = \^b) for spins on B sublattice. The phase 
diagram is shown in Fig. m 

At positive values of J z , the phases can be character- 
ized by the staggered magnetic and superfluid order pa- 
rameters: mf = \\S\ — 5fj|, and = ^l^i" — S^j, re- 
spectively, where Sj- = (S^, S^). The partially and fully 
polarized axial (or Ising like) antiferromagnetic phases 
(Al and A3, respectively) exhibit finite staggered mag- 
netization. Additionally, superfluid phases appear that 
break the U(l) (or 0(2)) symmetry: In the planar super- 
fluid phase SFq (i.e. the easy-plane Neel ordered phase 
realized in the Ba2CoGe207) only the fluid order param- 
eter 0[/(i) is non-zero, while in the conical superfluid 
phase SFa both the mf and 0[/(i) have finite expecta- 
tion values. In the conical ferromagnetic phase the finite 
magnetization along the z direction coexist with a finite 
Oj/(i) . A first order transition line separates the planar 
superfluid phase SF from the A3 and F3 gapped phases 
for smaller value of A, while for larger values of the on-site 
anisotropy, the SFq is bordered with the canted super- 
fluid states via a second order transition line. 

We remark that the J z -B- — J z symmetry of the phase 
diagram can be understood as follows: assuming copla- 
nar spin structure in a plane S perpendicular to the 
easy plane, rotating the spins on one of the sublatticc 
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FIG. 1: (color online) Variational phase diagram for h = as 
the function of A/ J and J z /J- Solid lines stand for continuous 
(second order) phase boundaries, while the dashed lines de- 
note the first order phase boundaries. The black point repre- 
sents the SU(2) symmetric isotropic Heisenberg limit. The Al 
and A3 label Ising-like antiferromagnetic phases, with spins 
S z — ±1/2 and S z = ±3/2 aligned along the z axis, while 
the F3 is the Ising-like ferromagnetic phase. The SFo is the 
easy plane Neel phase, the SFa is a conical Neel phase, and 
SFf is a conical ferromagnetic phase. 



by an angle tt around the axis perpendicular to the E, 
we change the sign of the J z , the conical antiferromag- 
net SFa changes to a canted phase, while the J and A 
remain unchanged. Thus we can map the J z > part of 
the phase diagram to the J z < part, as seen in Fig. [T] 
Within the variational approach, the large A stabilizes 
the easy-plane Neel-state (i.e., SF planar superfluid) 
ground state even for relatively large exchange anisotropy 
(J z /J J$ 4), with the following trial wave-function: 

|*a> = e-^^sp), 

|* B > = e-***> SF >, (2) 

where the I^&sf) is characterized by the single variational 
parameter rj, 

,^ sf) = II) + ^1) + ^]- + 

v/6?? 2 ± 2 

The angles ipA and (ps measure the tilting of the spins 
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from the x-axis and can be written as 



<Pa = 2 



^ = ~2 



(4) 



so that ip a — tfB — 7T- This describes antiparallel spins 
(Neel-state) on the two sublattices, as the expectation 
value of the spin components on the A sublattice is 



377(77 + 1) 
3f7 2 



1 



(cos^,sinv3A,0) , (5) 



with an analogous expression for the spins on the B sub- 
lattice. For 77 = 1, the a) and \^b) are spin-coherent 
states, with spin length 3/2. As the parameter 77 in- 
creases, the supression of the | ±3/2) states reduces the 
length of the spin, and the wave function describes a 
spin of mixed dipolar, quadrupolar, and octupolar char- 
acter. The independent parameter <j) in Eq. Q carries the 
0(2) symmetry breaking property of the superfluid phase 
(note that H commutes with S z ) and it is the direction 
of the spins with respect to the y-axis. For convenience, 
we choose <j> = hereinafter. 

The expectation value of the energy per site in the 
easy-plane Neel-state, as the function of parameter 77, 
reads 



E SF »(r)) 3 J) 2 



N 



4 3?7 2 



-A 



1 



1 v (v + ir 

(3t7 2 + l) 2 



J 



(6) 



where N denotes the number of spins. The E SF ° (77) be- 
comes minimal when the condition 



A _ 3(tj 2 - 1)(3?7+1) 
J ~ 



3?7 2 



1 



(7) 



is fulfilled. For only the in-plane spin components S x 
and S v are finite, the energy of the planar superfluid 
phase (E SF ° (77)) does not depend on the J z term, which is 
reflected in the 77 being dependent on the A/ J only. The 
usual procedure is to solve for the value of 77 as a function 
of the parameters of the Hamiltonian — in our case it 
would mean finding the solution of a cubic polynomial. 
In the following we will rather use the parameter 77 in the 
expressions instead of A as given by Eq. (0 , and replace 
the relevant 77 that corresponds to a specific A value only 
at the end. 

The first order transition between the antiferromag- 
netic A3 or ferromagnetic F3 Ising phases and the easy 
plane-Neel antifcrromagnct happens when E SF ° is equal 
to E A3 ' F3 = 9A/4 =F 9J z /2 (upper sign corresponds to 
A3, lower to F3 Ising phase), conditions that provide 
the phase boundaries 



±J 



4t7 3 (3t7 2 + 277-I) 



(3r/ 2 



(8) 



Let us briefly comment on the two limiting cases of the 
single- ion anisotropy. (i) When A/ J — > the 77 — 1, and 
the solution is the spin-coherent state mentioned earlier. 
Furthermore, for J = J z the Hamiltonian ([1]) recovers 
the 0(3) symmetry and the ground state is the spin-3/2 
Neel-state. (ii) In the limit A/J — > 00 the variational 
parameter 77 — > 00, the S z = ±3/2 states of the spins are 
suppressed and the ground state wave function is com- 
posed of the S z = ±1/2 states. Allowing for a general 
wave function in this reduced (two dimensional per site) 
Hilbcrt space, the tip of the spin spans a surface of an 
oblate ellipsoid. The length of the spin is maximal (equal 
to 1) when it lays in the xy-plane and minimal (equal to 
1/2) along the z-axis, therefore a finite antifcrromagnetic 
exchange selects ordering in the xy-plane. 



III. MULTIBOSON SPIN- WAVE SPECTRUM 

The usual spin- wave theory is a 1/S expansion in the 
length of the spin, where a single Holstcin-Primakoff bo- 
son is introduced to describe the transversal fluctuations 
about the classical ground state. As we shall see shortly, 
the multiboson spin-wave approach supports the inclu- 
sion of more bosons that have essentially different nature, 
providing a powerful method to discuss higher order, i.e. 
quadrupolc- or octupole-type excitations^"^ 

Let us first introduce the bosons aj^ that create the 
S z = m states of the S = 3/2 spin: \m) = a } n \ vacuum). 
Using the four a bosons, the diagonal spin operators can 
be written as 

3/2 3/2 

S Z = ma li a m, (S") 2 = m2a m a m' 

m=-3/2 m=-3/2 



and the off-diagonal spin raising operator is 



(9) 



s+ = V3 (4 



3/2 a i/2 



a -l/2 a -3/2 



)+2a\ /2 a_ 1/2 , (10) 



while the spin lowering operator can be easily obtained by 
its hermitian conjugate. All the spin operators (including 
higher order polynomials) can be expressed as quadratic 
forms, so that they keep the number of bosons M on 
each site conserved (M = J2 m a ln j a m ji an d M = 1 for 
the S = 3/2 spin). As a consequence, the Hamiltonian 
([T]) also commutes with the J2 m a m j a m j- Furthermore, 
written in this form, all the operators obey the expected 
commutation relations. 

For we want to carry out the multiboson spin- 
wave method starting from the planar antifcrromagnetic 
ground state, we apply an SU(4) rotation in the space of 
at, bosons: 



4 



1 



9 = 

J - 



v/6r/ 2 + 2 
1 



V lV _ 8 , + 2 L ^ ( e ^-3/ 2 " ^^/J + - 1) (e*W_ 1/2 
1 



1/2 



^/14t7 2 — 8?7 + 2 



(2^7-1) ( 



/2 c "1/2 



(11) 

(12) 
(13) 
(14) 



In this rotated basis the variational ground state given 
by Eq. ((2]) corresponds to the \^>a) = a] | vacuum) with 
<Pj = Va = tt/2 for spins on A sublattice and \^b) = 
a^|vacuum) with <pj = <ps = — 7v/2. The b, c, and d 
are suitably chosen bosons that will play the role of the 
Holstcin-Primakoff bosons and describe the excitations 
of the system. Namely, inverting Eqs. (fTl"T) -(fl"4" |) we can 
express the spin operators [Eqs. © and (|TU|)] using the 
a, b 7 c and d bosons, and replacing 

«J \I M b > 3 ^ ( 15 ) 

°j V 17 - b lh c h ^ ( 16 ) 

the spin operators (see Appendix [5} still follow the ex- 
pected commutation relations, analogously to the famil- 
iar Holstein-Primakoff transformation that uses a single 
boson. Performing an expansion in the parameter 1/M, 



one can further follow the procedure of the conventional 
spin wave theory for ordered magnets. The multiboson 
spin-wave Hamiltonian up to quadratic order in bosons 
then reads 

H = M 2 H {0) + M 3/2 H {1) + MT-L (2) + 0(M 1/2 ) , (17) 

where is equal to the variational (or, equivalently, 
the mean field) energy Eq. (|6]) and the is identically 
zero when the variational condition ([7|) is satisfied. The 
TiP^ is quadratic in bosonic operators and can be written 
as a sum of Hamiltonians in the k space, 

keBZ 

where the wgk = W$-k and <k = <l k read 



H 



(2) 
bd.k 



H 



(2) 
c,k 



6(77 + l) 2 (9?7 3 - 5?/ 2 -r? + l) 3 (7t7 2 - 4r/ + l) 

J — : ^7k + 



12?7 2 (?7+ If 



(3t7 2 + 1) (7?7 2 - 4?; + 1 
12t7 2 ( 77 + 1) 2 



+ 



(377 2 + 1) (7t7 2 - 4?7 + 1) 
72r7 3 (?7 + l) 2 



3tj 2 + 1 

3 (7T7 2 -477 + I) 
2 + 37/2 + 1 



(3?7 2 + 1) (7r7 2 -477 + I) 



7k(^L k + 6 k 6_ k 



(3t7 2 + 1) (7t7 2 - 477+I) 



j ( 4 d k + 



k"-k 



-6V3 



77(77 + 1) (77 -l) 2 
(3/7 2 + 1) (7i7 2 — 477 + 1) 
9(77- l) 4 



(677J - j, 7k ) ( &]a + &i k d_ k + 4&k + rft -k&-k ) . 



(3t7 2 + 1) (7t7 2 - 477+I) 
6^377(77 + 1) (77 -l) 2 



JzJw did 7 . + rf k d_ k + did. + cf k c2 



k -k 



(37? 2 + 1) (7t7 2 - 477+I) 



T) J,7k (& k d-k + &_,A + 6 k dL k + bl k d 



,t 



6J(77 + 1) (c k c k + c t _ k c_ k ) - ^te-lM^ 1)2 J 7k ( c t c L k 



&k\ + &V-k, 



(19) 
(20) 



In the above Hamiltonian we replaced A by the expres- sion for the 77, Eq. ([7]). The 

bi = ^=y, ^b[. 

b k = -LVe^ft, , (21) 
k v/A 4^ J 



5 



is the and b bosonic operator in the momentum space, 
with analogous equations for the c and d bosons. Setting 
the lattice constant to 1, the geometrical factor 7k can 
be expressed as 



1 



7k 



(cos k x + cos k y ) 



(22) 



Note that 7k = — 7Q+k ; where Q = (ir,ir) is the Neel 
ordering vector. Furthermore, for k — > 0, 2yT — 7k — > 
|k| = k, similarly 2^/1 + 7Q+k -> |Q + k| as k — > Q. 

The spin-wave Hamiltonians (flT)|) and (|20p can be di- 
agonalized using Bogoliubov transformation. The excita- 
tion spectra obtained by numerical Bogoliubov transfor- 
mation are shown in Fig. [5] for a set of selected A and J z 
values. For a given value of 7k we get three eigenvalues 
(i.e., six modes for each k in the reduced Brillouin zone, 

i 2) 

as shown in the figures), one from H. c (we will denote 

(2) 

this mode by letter 'c') and two from T-L bd (the 'b' and 'd' 
mode). This notation for the modes can also be traced 
back to the A — > limit, where these modes stem from 
the r] -+ 1 form of the bosons b, c, and d [Eqs. ([T2" |) -([n ) ] 
(see Appendix iBl for the discussion of the A-)0 limit). 
The 'b' band is the lowest in energy and goes linearly to 
at the k = (0,0) wave vector in the reduced Brillouin 
zone. The 'c' and 'd' bands are both gapped and higher 
in energy than the 'b' band, typically 'd' being the high- 
est. Their dispersion is much smaller than that of the 
'b' band, and disappears as we decrease the single ion 
anisotropy A, becoming flat (localized) for A = 0. 

(2) 

The eigenvalue of the "He can be calculated analyti- 
cally. It is independent of J z and reads 



6J(rj 




fa -1)2(37?+ I) 2 
1 : : ~7k 



(77 + 1) (3772 



1) 



(23) 



(2) 

The analytical expression for the eigenvalues of n^J is 
beyond our reach, except for two special cases: (i) along 
the lines k x + k y = tt, when 7k becomes zero and the 
energies are: 



Ub,d 



3J(?7 + 1) 2 (377 + 1) 



3?y 2 



1 



^ 3 J (7] + l)y/977 4 - 24t7 3 + 22t7 2 + 877 +T 
± z o (24) 



3 V 2 



1 



(ii) at k = (0,0) the 7k = 1, and one of the cigenmodes 
is the uib — Goldstone mode associated with turning 
the order parameter in the xy plane and desribed by the 
self-adjoint operator bU x „ = b r 



6 t xb 1 +b V3(ri-l) 2 U 
GM (0,0) (0,0) 277(77 + 1) V (°.°) 



>,o) 



, (25) 



that commutes with the spin- wave Hamiltonian"H (2) . At 
the same time, the 'd' branch has energy 

= IS^V^^TWl) ■ (26) 
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FIG. 2: (color online) Multiboson spin-wave dispersions 
in the reduced Brillouin-zone for (a) different single-ions 
anisotropies at J = J z and (b) different exchange anisotropics 
when A = 8J. The labels 'b' and 'd' denote the two branches 
that are eig envalues of the W^ ) k [Eq. (15)]. The 'c' (dashed 

lines) labels the eigenvalue of the H^l. [Eq. (20) ] which is 
independent of the value of the exchange anisotropy J z . 



In both of the above cases the energies are independent 
of the exchange anisotropy J z . Apart from these special 
instances, analytical result are available in the A — > 
and A — > +00 limits, as shown in Apcndices [Bl and [U] 

In Fig. [^b) we show the evolution of the dispersion 
as we change J z , while we keep A constant (A = 8 J in 
the figure). It is the 'b' band that is the most sensitive 
to value of J z , while the higher energy 'd' is only weakly 
affected. The energy of the 'b' mode is linear in mo- 
mentum, loi, = Vbk as k — > 0, where the velocity can be 
calculated analytically by expanding the Hamiltonian in 
y/1 — 7k and reads 

6r?(r7 + l) r 



Vb 



4?7 3 



3rf 



1 



77 3 — ?7 2 + 377 



1 



J + Jz 



(27) 



Increasing J Zi the energy of the 'b' mode decreases at 
7k = — 1, and becomes at a critical value J z = J%, 
where 

4„3 

J z = 5-^; 7- (28) 

ti 3 - if + 377 + 1 y ' 

The behavior of the 'b' mode can be connected to the 
phase boundary of the easy-plane Neel phase: its soften- 
ing marks the line of the second order transition into the 
conical Neel phase (or superfluid SF^), as shown in the 
variational phase diagram, Fig.[T] Quite interestingly, for 
ferromagnetic J z Ising coupling, the second order phase 
transition line into the conical canted phase that is given 
by J z — —J z , is indicated by the vanishing of the spin 
wave velocity Vb, Eq. (|2"7|) . 
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In order to test the reliability of the multiboson spin- 
wave method, we calculated the expectation number of 
bosons in the ground state, (b^b + c^c + d)d ). Wc 
found that the quantum fluctuations are the greatest and 
the boson expectation value are the largest in the fully 
isotropic (A = 0, J z = J) case: (b% + e'e + <vd ) = 
0.197. Introducing even a small anisotropy reduced this 
value considerably Our result for the isotropic point 
coincides with the known result for the spin reduction 
AS* = 0.197 in the square lattice,— as the bosons c and 
d decouple from the system, and the boson b plays the 
role of the standard Holstein-Primakoff magnon. 
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IV. SPIN STRUCTURE FACTOR 

The intensity at energy u and momentum k in inelas- 
tic neutron measurement is determined by the magnetic 
cross section 



rf 2 cr(k, uj) 

OC 7 

dujdTL ^ 



Snv — 



(29) 



at zero temperature (fi, v = x 1 y,oi z), where 



S">,k) = 52(0\S» k \f)(f\SZ\0)6(u> - Uf ) (30) 
/ 

is the dynamical spin structure factor, |0) is the ground 
state, and the summation is over the excited states /. 
The matrix elements and energies in the S' A ' l/ (u;,k) can 
be calculated using the spin wave theory. The 1/M ex- 
pansion for the spin operators, using Eqs. ((9]), (fTUf and 
Eqs. dHJ-dHl), reads 



±M 



±VM 



i\/3 y/77] 2 -477+I 
2y/3 v 2 + 1 
3r?(77 + 1) 
3t7 2 + 1 
V3(r? 



M, 



l)(3r? + l) 



2 (3r? 2 + 1) 
V3?7(?7 



(31) 



(32) 



v/377 2 + lv/7?7 2 - 477 +l 



3(r7 -l) 2 



2y/3rf + li/V — 477 + 1 



d] +dj 



(33) 



where only the leading order terms that are proportional 
to M and y/M are shown, and the upper (lower) sign 
corresponds to the spin on sublattice A (B). Due to the 
alternating sign in and Sj , in the Fourier transform of 
these operators the Holstein-Primakoff bosons are shifted 
by the Neel-ordering vector Q = (tt, 7r): 



SI cx -i(^+Q _6 -k-Q 



M 



k+Q 



-k+Q 



M 



(34) 
(35) 
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FIG. 3: (color online) Dynamical structure factor S' MM (k, uS) 
with /i G {x,y,z} for J z /J = 1 and (a) A = J/2 and (b) 
A = 6 J along a path in the Brillouin zone. The widths of the 
filled curves above the excitation energies (black solid lines) 
denote the strength of the matrix elements |(/|>S'k M |0)| . The 
S xx (k,iu) diverges as 1/lu as k — >■ Q. We note that the ground 
state is a Neel antiferromagnet with the spins chosen to be 
parallel to the y-axis. Therefore, the 'c' is a stretching (or 
longitudinal) mode, associated with length fluctuations of the 
spins, while the 'b' and 'd' are transverse modes. 



correspondingly 
S""(w,k) = 



^|(/|5^|0)| 2 %-a; k+Q ) 



(36) 



for /j = x, y. Such a momentum shift is not needed for 
the S£. 

The matrix elements of the different modes in the Bril- 
louin zone of the square lattice, evaluated numerically, 
are shown in Fig. [3] for a small and a large value of the 
single-ion anisotropy A. It turns out that the 'b' and 
'd' modes have finite matrix elements with the S x and 
S z spin components that are perpendicular to the ori- 
entation of the spins, i.e. these modes are transversal 
modes, similarly to the modes in the conventional spin- 
wave theory. The 'c' mode is more interesting, since the 
only nonvanishing matrix element is with the S v spin op- 
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(ji/2,ji/2) (0,0) (ji,0)(ji/2,ji/2) (0,0) (ji,0)(ji/2,ji/2) (0,0) (lt,0) 



(k x ,ky) (k x ,k y ) (k x ,ky) 

FIG. 4: (color online) The different components of the dy- 
namical structure factor are shown separately in the reduced 
Brillouin zone for the anisotropy parameter A = J/2 in (a)-(c) 
and A = 8J in (d)-(f), while in all cases J 2 /J = 1. Here the 
widths of the filled curves above the excitation energies denote 
the matrix elements | {/|S'^ M |0) | 2 multiplied with energy ujf, 
i.e., oj5'' 1M (k, u}). The ground state is a Neel antiferromagnet 
with the spins chosen to be parallel to the y-axis. The dotted 
lines denote silent modes in the corresponding S MM (k,w). 



erator: in this mode the length of the spin changes ("spin 
stretching mode"). In all cases the S xx (k,uj) diverges as 
1/lo when k approaches the Q = (7r,7r) ordering wave 
vector — this reflects the zero energy cost of rotating the 
spins in the easy plane, i.e. the Goldstone-mode. On the 
other hand, the S zz (k, oj) associated with spin fluctua- 
tions perpendicular to the easy plane, while finite at Q 
for finite values of A, diverges as the anisotropy gap is 
closed. Eventually, S xx and S zz become equal for A = 0, 
when the full 0(3) symmetry is recovered. Furthermore, 
as A is decreased, the intensity of the 'c' and 'd' modes 
decreases rapidly (see Appendix IB"! for a detailed discus- 
sion of the A — > limit). To eliminate the 1/uj divergence 
and obtain a better comparison of the matrix elements, 
we show uiS^^Ck, uj) in Fig. (j4]), this time along a path 
in the reduced Brillouin zone. 



V. COMPARISON WITH Ba 2 CoGe 2 7 
NEUTRON SCATTERING EXPERIMENTS 

Inelastic neutron scattering measurements on 
Ba2CoGe207 were reported in Ref. [3 In the ex- 
periment constant-momentum scans were performed 
up to 3-4 meV (depending on momentum), and the 
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h (r.l.u.) 

FIG. 5: (color online) The spin-wave dispersion for A = 1.15 
meV, J z = 0.154 meV, and J = 0.197 meV. The wave vectors 
q used in Ref. [HI and in this figure for easier comparison are 
defined using the two-Co ion unit cell, so that q = (k x + 

ky , k x ky , 0)/27T. 



peaks in the intensity were traced to get the dispersion. 
A single mode was observed with a dispersion of 2.2 
meV and with a large anisotropy gap comparable 
to the dispersion itself. The mode was fitted using 
'conventional' spin-wave theory based on the Heisenberg 
model extended with strong exchange anisotropy only 
(i.e. A = 0). No higher modes were observed in the 
aforementioned energy window. 

On the other hand, additional peaks have been ob- 
served in far infrared spectral at around ITHz beside 
the 0.5 THz mode that corresponds to the 2.2 meV peak 
detected by inelastic neutron scattering (we use that 
IThz « 4.13meV). Furthermore, a recent study of the 
far infrared absorption in high magnetic field with high 
resolution revealed several modes that could be described 
by the Hamiltonian ([1]) using the presented multi boson 
models. Therefore we anticipate that those lines shall 
also be present in the inelastic neutron spectra. To be 
more precise, in Fig. [5] we compare the calculated multi- 
boson dispersions using the parameters A = 1.15 meV, 
J z = 0.154 meV, and J = 0.197 meV of Ref. ]M with the 
inelastic neutron scattering peaks taken from Ref. [l4|. As 
shown, the k dependence is nicely reproduced, and we ex- 
pect the additional, weakly dispersing peaks with smaller 
intensity at energies that are about 4 meV. 



VI. CONCLUSIONS 

In this paper we have extended the multiboson spin- 
wave theory to the spin-3/2 Neel antiferromagnet with 
strong single-ion anisotropy: the spin operators in the 
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ordered state are described by three Holstein-Primakoff 
bosons, and accordingly, the excitation spectrum in this 
approach consist of three modes for each spin. 

In the absence of the single-ion anisotropy, only one 
mode out of the three has finite matrix elements with 
the spin-dipole operators. This mode is equivalent to the 
magnon mode of the 'conventional' spin-wave theory in 
every aspect. The additional two modes, silent in spec- 
troscopy probes that interact with the magnetic moment 
only (such as neutron scattering), describe quadrupolar 
and octupolar fluctuations of the S — 3/2 spin and may 
become visible when these multipolar fluctuations couple, 
for example, to electric polarization. This may happen 
if the spins are in a non-centrosymmetric environment, 
where coupling between the quadrupolar spin fluctuation 
and electric polarization makes them visible in light ab- 
sorption experiments, such as in far infrared spectr a 19 ' 20 . 

The picture above changes markedly when finite single- 
ion anisotropy is present in the system. As a result of the 
single-ion anisotropy, the spins in the mean field (equiva- 
lently site-factorized variational wave function) approx- 
imation are not any more spin-coherent states in the 
Neel ordered phase: the suppression of the S z — ±3/2 
spin-states for A > leads to shortening of the spins. 
This allows longitudinal fluctuations of the spins (stretch- 
ing modes) that have finite matrix elements with the 
spin-dipole operator parallel to the ordered moment, 
thus become observable in neutron, electron spin reso- 
nance and other spectroscopics. Our findings are closely 
related to the observed longitudinal excitation mode 
in the pressure-induced in-plane AFM phase of N1CI2- 
4SC(NH2)2^ A similar approach has been introduced 
to describe the nature of the excitations in the gapped 
TICUCI3 spin-dimer compound, where the one-magnon 



Raman scattering has been found efficient for selectively 
observing such longitudinal excitations in the pressure- 
induced ordered phasesj 25 ' 26 

Finally, we compared the multiboson spin- wave modes 
with the inelastic neutron spectra of Ba2CoGe2C>7. We 
found that the lowest mode reproduces the measured dis- 
persion using the exchange and on-site anisotropy param- 
eters fitted from the evaluation of the far-infrared absorp- 
tion measurements^ - We propose that the other, higher 
energy modes shall also be observed in neutron scatter- 
ing experiments at «4mcV, with weak dispersions. From 
the analysis of the spin-structure factor the stretching 
modes could, in principle, be identified. Furthermore, we 
believe that such stretching modes shall appear in any 
S > 1/2 material with strong single-ion anisotropy, start- 
ing from the related compounds Ca a; Sr2- x CoSi20 | 2?l2S 
and Ba 2 MnGe 2 7 2°.. 
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Appendix A: Expansion of the spin operators 

Considering the bj, Cj, and dj as Holstein-Primakoff 
bosons and replacing the a' and a bosons with the ex- 
pressions (fT5|) and (Til))) , the spin operators on the A sub- 
lattice read 



iV3\/777 2 -477 + I 



i (3?7 2 



M 



2r) - 1) 



v/3?7 2 + 1^7t] 2 - 4?7 + 1 



b ?3 



b\c- 



c j c j 



M - bib, - etc, - d\d. 



33 



2VV- 4? 7 



t(«5< 



33 



d'e 

U 3 C 3 



3 3 



(Al) 



M 



3r/(r; + 1) 
" 3t7 2 + 1 
V3 (3?7 2 -277-I 



+ 



(2T7 2 



6?7 2 
+ 77 



7?7 — 4?7 



M - 

73 (577 



( -C ■ 



d 3 d 3 C 3 



677 



Uif -877 + 2 



Hi* 



^3 



M 



cc ■ 

3 3 



¥3 



3(1 - 377)77 x 



(1 ~ 37?) ctc.+ did- 
3 ? 7 2 + l ^ 7t7 2 - 477 + I J 



(A2) 
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s * = ., -T ,. , ( ^ M b h c h 4 d A + d W M V', c Is - d 'A 



2^3r/ 2 + l^V - 477 + I 

V^ + l) f./Vr X 7 



^ V + 2 " - ^ + 40 - .. ,-^, n . + . (A3) 



and 



2^3772 + VV -477+I V J J J Jy v/37J 2 + ly/7r} 2 — 477+I 



12?7 2 + 4 V ■? J 3 3 3 3 



3t] 2 + 1 



M h h c h <>:<>/:, + sV M - 6 & - c h d ) d 3 



(31772-477 + 1) ^ + 2^77(277 - 1) , + + (2V + 1), + (3 V - 88, + 9) 



28t7 2 - I677 + 4 J J 7t/ 2 — 4?7 + l V J J J J / 12?7 2 + 4 ' J J 28t7 2 — I677 + 4 33 

As written, the expression above satisfy the exact commutation relation of the spin operators, i.e. no approximation 
has been made yet. 

The spin operators expanded in 1/M, neglecting the 0(M -1 / 2 ) terms, are: 



2y/3r] 2 + 1 



* + ^ *) ( b \ c . _ c t^ - ^ 



5 / 



^7 



VVTTVV- 477 + 1 ^ ^ ' 2^2-477 + 1 ' (A5) 

3( ";"' u + 4) - ^ ( bj + tj) 

2^377 2 + V777 2 - 477 + I V 3 3 J ^3772 + 1^7772-477+1 V J V 



2^/3772 + 1^/7772-477 + 1 V J J J V ^3t7 2 + V777 2 - 477+I V 3 3 
(sz? _ 3(77 2 + 3) 2^377 ~, a 3 (772 + 3) , 

J } " 12772 + 4 M + 3772 + 1 V ^ V 5 ' + C ' 3 ) 1277 2 + 4 Vi b i + C3 ° 3 + d i d i 

(31,2-477+1) t 2^377(277 - 1) / t t s + , (39^- 3677 + 9) 

28772 - I677 + 4 3 3 + 7772-477 + I + a A 3 ) + 12772 + 4 + 28772 - I677 + 4 *vV 1 j 



I 

Appendix B: The case of small single-ion anisotropy Substituting this expression into the multiboson spin- 
wave Hamiltonian (|B2|) . we get 0(A 2 ) 



When A <C J, Jz , the parameter 77 in Eq. (0 that min- 
imizes the energy can be expanded as +3 (J + J z ) (^k^-k + b k b 



«$ k = I 6 J - 3 (■/ - Jz) 7k] + 

2-7k-^7k) (^k + ^-k) 



A 
2" 

a r j z 



. l-y ]7kK&L k + & k &_ k 



? ' = 1 + 67 + TSj5 + ° (jJ ) • (B1) + 18 J (4^k + dl k d_ k ) , (B2) 

«S = ( 1 2^ + A)fc k c k + cL k c_ k N ) . (B3) 
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The &, c, and d bosons decouple in this limit. The b 
bosons propagate, while the c and d bosons are localized. 
The energy of the different modes (here we keep higher 
order terms in Eq. (|B2[i ^ are 



6 Jy/1 - 7k 



•^ll + i^ + m + iJ-Jzi*)-, 



J 3J 



12J + A+(l-2 7k ) 



72 J 3 



A 2 
24J 



18 J 



A 3 
48J 2 



(B4) 
(B5) 

(B6) 



The uib + shows the typical square-root behaviour of the 
anisotropy gap on the exchange anisotropy J — J z and 
single-ion anisotropy A. The c and d bosons are related to 
higher, AS y = ±2 and AS y = ±3 transitions in systems 
(quadrupolar and octupolar), as for ?y = 1 the 



2c h 



■3d& 



(B7) 



For A = there are only bilinear exchange terms in the 
Hamiltonian, therefore the quadrupolar and octupolar 
moments of the spin on neighboring sites do not inter- 
act. The energy of the AS y = ±2 'c' mode is simply 12 J, 
and that is actually the Zeeman energy that corresponds 
to this transition in the mean field of the four neighbor- 
ing antialligncd spin, each contributing by (3/2) J to the 
Weiss field /iweiss = 4x3 J/2 = 6 J. Correspondingly, the 
energy of the <i-branch is 18 J = 3ft.\veiss, the Zeeman en- 
egy of the AS V = ±3 transition. Both 'c' and 'd' modes 
acquire dispersion in higher order in A: the mode 'c' in 
the order A 2 / J and 'd' in the order A 4 / J 3 . 

The matrix elements in the dynamical spin structure 
factor are given in this limit as 



= \\j J jtj + lk m +0{k2) (B8) 



l(o|^)l 5 



l^|c)| 2 




J - Jjk 



48J 2 

while both \(0\S*\d)\ 2 and \(0\S 



Jzlk + A/3 
0(A 3 ) 



0(A 2 ) (B9) 



(BIO) 



arc <3(A 4 ). The 



|(0|S'^|6)| 2 matrix element of the 'b' mode, associated 
with the Goldstone mode, diverges as 1 /w k +Q as we ap- 
proach the ordering vector Q. Similarly, the |(0|5^|6)| 2 
is also oc l/wk, however the divergency is now cut off 
by the anisotropy gap. In the absence of anisotropy, 
S xx (k,uj) = S zz (k,uj), as expected. As for higher en- 
ergy features, the weights of the 'c' and 'd' modes in the 
spin response vanish with A 2 / J 2 and A 4 / J 4 , and become 
negligible small for small values of A, as we already no- 
ticed in Fig.[3£a)-(c). 



It is also instructive to compare with the conventional 
spin- wave theory, where we keep the b boson only (setting 
r] = 1 in Eq. ([T2|) is exactly the form of the spin operators 
using b as the Holstein-Primakoff boson in the conven- 
tional case). The spin- wave Hamiltonian then reads 



H {2) = S[4J-2(J-J z ) lk ^ 
+S[2(J+J z ) lk -A] 



A] 



&3a 



-k u -k 



^-k 



M>-k (BH) 



Replacing the 5 = 3/2 value, the terms proportional to 
J and J z are identical to the ones in Eq. (|B2|) , while the 
coefficients of the terms involving A differ. 



Appendix C: Modes for large single-ion anisotropy 

From Eq. ([7]) we can obtain the 1/ A expansion of r/ in 
the A > J, J z limit: 



V 



A 
3J 



4J 

T 



o 



A 2 



(CI) 



Diagonalizing the Hamiltonians p^|) and (fT9"| in this 
limit, we get the following excitation energies 



4 = V(l - 7k) J (4 J + J 27 k + 0(1/A)), (C2) 



w k = 2A + 4J-3J 7k + 0(l/A), 
w£ = 2A + 4J-3J 7k + 0(l/A), 

and the matrix elements are 



|(0|^|6)| 2 



i<omi : 



Vij- 


Jzlk 


2^74 


- Jlk 


y/J- 


- Jlk 



2^4 J + J zlk 



0(1/A) 
fO(l/A) 



|(0|^|d)| 2 
\(0\Sl\d)\ 2 



qr . 3_^, +0(1/A2) 
3__^^ +0(1/A2) 



3J 2 
A 2 



(l-7^) 2 + 0(l/A 3 ) 



(C3) 
(C4) 



(C5) 

(C6) 

(C7) 
(C8) 
(C9) 



The weakly dispersing 'c' and 'd' modes describe excita- 
tions of the | ± 3/2) S z states. We see that almost all the 
matrix elements arc finite when the single-ion anisotropy 
is large. The smallnes of the |(0|S£|d)| 2 can be under- 
stood as these are the fluctuations that need to overcome 
the large A energy. 

Here we can seek the correspondence to an effective 
model. When A is large, the S z — ±3/2 states can be 
projected out from the Hamiltonian. The essential de- 
grees of freedom are reduced to the S z = ±1/2 states of 
the S = 3/2, and these two states and the interaction be- 
tween them can be represented by the a z = ±1/2 states 
of an effective spin-1/2 model with the following XXZ 
Hamiltonian: 

n cfi = ( 4J K<\" + *M) + J ^M) ( C1 °) 

<i,3> 
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where the of are spin- 1/2 operators on site i that act on netic spin- 1/2 problem that provides the same dispersion 
the Hilbert space of the effective spins. So in the leading as in Eq. (|C2[) . 
order we are left with the usual anisotropic antiferromag- 
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